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In lattice QCD spectrum calculations, it is desirable to obtain multiple excited state energies in 
each symmetry channel. Typically, one constructs several interpolating operators for the sym- 
metry channel of interest, forms the 'correlator matrix' of all possible two-point functions, and 
uses the variational method to obtain as many energy levels as possible. We present a detailed 
look at this last step, starting from a discussion of the symmetry properties and spectral struc- 
ture of the correlator matrix. We continue by motivating and describing the variational method, 
before discussing some conceptual and practical challenges concerning the light baryon sector. 
We conclude by mentioning some alternate spectrum extraction methods currently under study. 
Throughout, we attempt to quantify all approximations and assumptions, and we illustrate our 
points using a nucleon correlator matrix estimated on dynamical two-flavor lattice data. 
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1. Formalism 

We will work with a D x D matrix of two-point functions, the so-called correlator matrix: 



-(P-T)HMP) e -TH-p?(P) 



(1.1) 



where Z = e~$ H is the partition function, j8 is the temporal extent of the lattice, H is the Hamilto- 

ip\ 

nian of the system, and the D operators {€? a } are designed to have: 

• Zero three-momentum: p = 0, 

• Gauge invariance: singlet color structure, 

• Definite lattice spin [jl]]: G\,H,or G2 (double- valued irreducible representations 1 of Oh), 

• Definite parity P = ±, also denoted as: g = gerade (even) and u = ungerade (odd). 

The operators carry baryon number B = 1 and have been constructed to have the following behavior 
under charge conjugation, represented by the unitary operator Uc'- 
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where we have dropped the dagger because {Uc) 2 = I- This behavior leads to the following relation 
between the positive- and negative-parity correlator matrices: 
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where we have used the cyclic property of the trace and the fact that charge conjugation is a sym- 
metry of the theory: [H,Uc] = 0. This relation is used to improve our statistics by averaging our 

(p) (-P) 
numerical estimates of C ah (t) with those of C ha (j8 — t). 

The system is defined in a finite volume and thus has a discrete spectrum. Writing the trace as 

a sum over energy eigenstates {\n)} and inserting a second complete set of eigenstates {|&)} gives 

us the general spectral representation of the correlator matrix: 
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(1.9) 
(1.10) 

(1.11) 
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The operators we use do not mix among the different lattice spin irreps. 
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In Eq. II. Ill we may freely divide a factor of e °P out of the bottom and top, effectively cancel- 
ing any zero-point energy of the theory: H\0) = 0. Each of our baryon operators 6 ' a creates a 
three-quark system of a given parity P and annihilates a three-antiquark system of the same parity 
P. Unlike bosons, a fermion and its antifermion have opposite intrinsic parities; an antibaryon's 
energy is therefore the same as that of its corresponding opposite-parity baryon state, its so-called 
opposite-parity partner. Since chiral symmetry is broken in QCD, the masses of the opposite-parity 
partners may differ. 

(p\ (p\ 

Assuming throughout that (k\0 b \n) < (k'\€> ' b \n') for n,k < n',k', the leading terms in the 
spectral representation involve the following non- vanishing matrix elements: 

rp\ tp\ 

• {k\& a |0) where \k) is a baryon state (B = 1) having parity P and energy , 

tp\ l—p) 

• (0\ff a \n) where \n) is an antibaryon state (B = — 1) having parity P and energy E,, , 

ip\ 

• (BB\ff a \B) where \B) has baryon number B = 1 and \BB) has baryon number B = 2, 

ip\ _ _ _ 

• \B\ff a \BB) where \B) has baiyon number B = — 1 and \BB) has baryon number B = —2. 

In the last two matrix elements above, a positive-parity operator will connect single- and multi- 
(anti)baryon states with the same parity, while a negative-parity operator will connect single and 
multi-(anti)baryon states with opposite parity. 

Treating the largest terms explicitly, and using charge conjugation to convert the antibaryon 
states \n) into baryon states \k') (the trace is invariant under charge conjugation) gives: 



+ (l|^|l,l)(l,l|^|l)e- £|T ^ £l(/3 - T) 

+ •••. (1.12) 



where E\ is the energy of the lightest baryon state |1) accessible by the action of our operators 
on the vacuum, regardless of parity. We have approximated the energy of the lightest interacting 
B = 2 state |1, 1) as E\, to be conservative. Thus, all higher-order terms are suppressed by a factor 
of <?- £ '0, which is <9(lCr 6 ) in our studies. At this point, we note that before neglecting terms, it is 
important to keep an eye on e~ Ek% vs. e~ E[ P . Dropping the higher-order terms gives us a spectral 
representation of the correlator matrix in terms of single-baryon states only: 
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i — i (p) 

Notice that this approximation maintains the previously established relation (Eq. 1.8): C h (t 



C ba (/$ — t). In the following, we will treat the approximation (Eq. 1.13) as exact and will re 



fer to the two terms in Eq. |1.15| as the 'forward' and 'backward' parts of the correlator matrix, 
respectively. For convenience, we define: 



3?w 



(1.16) 



which allows us to write: 



(1.17) 
(1.18) 



It is our ultimate goal to extract the spectrum E) ',Ejj for the various QCD hadron sectors. 



2. The variational method 

Ideally, we would like to work with a very large temporal extent j6 >> T such that: 

C%\^C^{x) = {^V m Wf\ (2-1) 

in which case the problem reduces to simply diagonalizing the operator e~ Ux in the subspace de- 
fined by the trial basis This can be accomplished with the variational method [Q, [3| Q]. 

We begin by refining the states in our trial basis by (formally) evolving them by some amount 
To/2 in imaginary time. This serves to 'relax' them onto the lower-lying energy states of interest: 

= e- H ^ 2 \&P), (2.2) 
= e - Ht °/ 2 ( p\O), (2.3) 
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k=\ 



(2.4) 



where we have used the fact that (0|^ a |0) = 0. For sufficiently small To, the trial basis is linearly 
independent, but is not necessarily orthogonal or normalized: 

(^ P) (to)|< ) (to)}=C^(t ). (2.5) 

We can use the variational method to diagonalize the operator e~ // ( T ~ T o) in our trial subspace by 
forming the linear combinations: 

D 

\&i) = Y.\^ P) ( T o))v ai , (2.6) 

a=\ 

where the complex numbers v ai - are the variational coefficients obtained by finding the extrema of: 

(0,.| g -g(*-*b)|0,.) _ VJC( P )(T)V; 

vJcW(T )v,-' 
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Setting the derivative of Eq. 2.7 with respect to the v ai - to zero yields a generalized eigenvalue 
problem whose solutions are the columns v, of the variational coefficient matrix V: 



c( p )(T)v,- = A ( -d p )(ToM. 



(2.8) 



The eigenvalues {A,} appearing in Eq. are known as principal correlators, and if all has gone 

_p(p) / , 

well, they should asymptote in x as e > ^ T ~ T «" while the v k maintain the orthogonality of the 
associated states {| ©,•)}. From the properties of the generalized eigenvalue problem, we can see 
that the metric of the vectors v, used to diagonalize O p '(t) at later times is fixed at the reference 
time To by: 



(& i \&j)=v^ p \xo)v j = 8 i 



(2.9) 



It should be noted that the v a i, h, and |©,-) all depend on both x and To- 

tria 

'|0> 



From the spectral expansion of the D trial basis states in Eq. |2.4| , we see that if we can choose 

~ e -E k T /2 {k ^(P) t 



To large enough such that 



is negligible for k> D, then Eq. M implies that: 



(k\®j) = £ (k\e- H ^ 2 ff { p\0)vaj « 8 kj , j,k < D, 



(2.10) 



a=\ 



and that 



(& k \e- H ^\& k ) = v\6 p \x)v k =X k * e- E f'^ (2.11) 
occurs almost immediately for x > Xq. If Tq is not sufficiently large, it can be shown that the 



asymptotic form (Eq. [2.1 1| ) is still reached, but more gradually. In this case, it is possible that the 



(P) 

signal will be lost in the noise before a reliable estimate of E, can be made. 



We would like to choose To to be large enough such that Eq. |2.10| holds to within statistical 
uncertainty, but are limited by three main factors. First, we must estimate the correlator matrix 
numerically; the signal-to-noise ratio decays exponentially, leading to a noisy estimate of the trial 



basis for large To- Second, as we increase To towards j8/2, the spectral terms in Eq. |1.13| become ex- 
ponentially small, and thus the rank of C^ p )(To) will drop below D. When this occurs, the condition 
number of the matrix will blow up because the trial basis will not have full rank, and the generalized 
eigenvalue problem will be numerically unstable. This could be remedied by selecting some subset 
of operators (reducing D), but would incur the cost of throwing away information-containing en- 
tries of the original correlator matrix. Finally, as we increase To, the previously neglected backward 
term C^ p ^*(j8 — To) will begin to contribute, with the particularly serious consequence that: 



i P) (ToM P \<to)}^C%>(xo). 



(2.12) 



This means, for example, that as we increase To to reduce the positive-parity contamination, we 
will introduce negative-parity contamination. These issues are not always troublesome in practice, 
but we have found that they are indeed present in our nucleon spectrum study. 
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3. Eigenvalue behavior 

It is illuminating to examine the eigenspectrum (and rank) of an actual lattice QCD correlator 
matrix (Fig. [I]) as a function of x before applying the variational method. Here we consider a corre- 
lator matrix defined using D = 16 previously designed [JXJ, ^, ||] nucleon operators (isospin 7 = 1/2, 
^3 = +1/2), and estimated from 336 LHPC 2 two-flavor anisotropic Wilson lattice configurations 
having the following parameters: 

• Lattice dimensions: 24 3 x 64, 

• Spatial lattice spacing: a s w 0.11 fm, 

• Anisotropy: a s /a z 3.0, 

• Pion mass: m % 400 MeV. 

To find the D = 16 eigenvalues Ay as functions of X, we solve C^ p \x)vj(x) = A/(t)v/(t) at 
each x. The largest four A/(t)'s, corresponding asymptotically to the lowest four energy levels, 
are plotted on the left in Fig. [I]. The most prominent feature of this plot is the apparent crossing 
and transient kinking of the eigenvalue flow with X. The blue curve in particular exhibits interest- 
ing behavior: it falls as exp[— E% x] for early x, then becomes contaminated with the backward 
ground state exp[— e\ '(j3 — t)] at intermediate X, and finally is forced (by the higher eigenval- 
ues) to asymptote as exp[— E^ '(/j — x)]. We may understand this qualitatively by the fact that 
the negative-parity operators excite the negative-parity energy levels differently than the positive- 
parity operators excite the positive-parity energy levels. In other words: the coefficients diagonal- 
izing C' p '(t) for x « j8/2 are largely independent of the coefficients diagonalizing C^ F \x) for 
X >> /J/2. Unfortunately, our variationally-determined negative-parity operators are completely 
determined by our variationally-determined positive-parity operators. 




i / \ x / a t / a T 

Figure 1: These three plots are associated with a 16 x 16 nucleon correlation matrix C^ p \x) in the G\ 
lattice spin channel. Left: the largest four eigenvalues of the positive parity matrix. We see crossing and 
transient behavior due to the interference of the forward and backward parts of the correlator matrix. Because 
our negative parity operators have a fixed relation to our positive parity operators, we do not have enough 
variational coefficients in this technique to separately isolate the positive- and negative-parity contributions 
when they are both present. Center and right: the largest four principal correlators for the positive- and 
negative-parity correlator matrices, using To = 4a T and To = 5a T , respectively. In the center plot we see 
suspected crossing behavior around T = 30a T and again around T = 45a T . In the right plot we see significant 
opposite-parity partner contamination (from the proton) starting around T = 10a T . 



2 Special thanks to Robert Edwards and Balint Joo from Jefferson Lab. 



6 



The Spectral Structure of Correlator Matrices 



Adam C. Lichtl 



To understand this behavior, we look at the spectral structure of A/(t): 

Xj(T)=v](z)d P \z)v J {z) = v;(T)C( /> )(T)v y (T)+v}(T)C(- P )*( i 8-T)v ; .(T), (3.1) 

= V ; t (T)C( P )(T)vy(T)+vJ(T)C(- P )( i 8-T)v}(T), (3.2) 

ee {®f \e- H *\®f } ) + (e { r^\e- H ^\® { r P) ), (3.3) 

where we have introduced the states \®j } ) created by linear combinations of our basis operators: 

< = I &i F ^t^-\r 0-4) 



The variationally-determined ; operators behave as follows under charge conjugation (c.f. Eq. |1.2j ) 



Uc® { pU c = £ UcdpU&a] = £ d P) v aj = ® { j P \ (3.5) 

0=1 0=1 

and thus we see that a brute-force application of the variational method to a correlator matrix 
having significant contributions from both forward and backward components will be plagued with 
contamination because we are using only one set of variational coefficients v a j to diagonalize two 
matrices O p )(t) and 0~ p )*(/3 — z), which separately have different eigenvectors. 

4. Outlook 

We are currently examining two alternate techniques for the extraction of the spectrum from 
a correlator matrix. First, we are considering the possibility of introducing more coefficients by 
taking linear combinations of the elements of (t) and C^ p ^ Q3 — t). This may allow us to better 
decompose (t) by breaking the charge conjugation relationship of our variationally determined 
operators. A second approach would be to attempt full matrix fits of increasing rank M to the 
reduced matrix Rij(T) = v]c{x)vj, where the v,'s are variationally determined at some fixed time 
T* and only M of the column vectors v, are used (i,j = 1, . . . ,M) . We expect continuing progress 
towards producing reliable determinations of the various excited baryon spectra predicted by QCD. 
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